Schwinger boson mean field theory of the Heisenberg 
Ferrimagnetic Spin Chain 

Congjun WuT+j Bin Chen*, Xi Dai*, Yue Yu*, and Zhao-Bin Su* 
* Department of Physics, Peking University 
* Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100080, China 



Abstract 

The Schwinger boson mean field theory is applied to the quantum ferri- 
magnetic Heisenberg chain. There is a ferrimagnetic long range order in the 
ground state. We observe two branches of the low lying excitation and calcu- 
late the spin reduction, the gap of the antiferromagnetic branch, and the spin 
fluctuation at T = OK. These results agree with the established numerical 
results quite well. At finite temperatures, the long range order is destroyed 
because of the disappearance of the Bose condensation. The thermodynamic 
observables, such as the free energy, magnetic susceptibility, specific heat, and 
the spin correlation at T > OK, are calculated. The T\uni has a minimum at 
intermediate temperatures and the spin correlation length behaves as T _1 at 
low temperatures. These qualitatively agree with the numerical results and 
the difference is small at low temperatures. 
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I. INTRODUCTION 



A variety of exotic physical phenomena in the low dimensional magnetic systems have 
been attracting much interest in recent years. In these systems, the physical pictures ob- 
tained from the classical approach are often greatly modified or even in contradiction with 
the result of the strong quantum fluctuation and topological effect. Haldane [|IJ conjectured 
that the one-dimensional integer-spin chain with the nearest neighbor coupling has an en- 
ergy gap in the spin excitation spectrum and the spin correlation decays exponentially with 
distance, whereas that of the half-odd-integer spin chain is gapless and the spin correlation 
decays algebraically with distance. 

It is very interesting to discuss the physical phenomena for the spin chain mixed by 
different kinds of spins. Recently, one-dimensional Heisenberg ferrimagnetic spin chain, 
which is made of two kinds of spins, S A = 1/2 and S B — 1, has been considered. This 
one-dimensional chain can be described by the Hamiltonian: 

N 

H = ■ S i+V , (1) 

2=1,77 

where the antiferromagnetic coupling energy J is set to equal 1. N is the number of unit 
cells, and rj is the index of the nearest neighbors. Brehemer || et al showed that the 
absolute ground state of this model has a ferrimagnetic long-range order and obtained the 
low-lying excitation, by using the spin-wave theory(SWT) and the Quantum Monte Carlo 
method [QMC], which is confirmed by Kolezhuk || et al with the matrix product approach. 
It is also consistent with Tian's @ rigorous theorem that the absolute ground state for one- 
dimensional antiferromagnetic Heisenberg model with unequal spins has both antiferromag- 
netic and ferromagnetic long-range orders. Furthermore, Yamamoto et al |J used the mod- 
ified spin wave theory (MSWT), the density matrix renormalization group method(DMRG) 
and QMC to calculate the thermodynamic observables. 

In this paper we study the ferrimagnetic spin chain by means of the Schwinger boson 
mean field theory(SBMFT). The theory has been applied successfully to the infinite Heisen- 
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berg chain [Q , presumably due to the neglect of topological excitations in the SBMFT. It can 
also be extended to the case of the magnetic order range |10[ by identifying the magnetic 
order with the Bose condensation of the Schwinger bosons. In the Heisenberg ferrimagnetic 
spin chain, we find that the SMBFT theory is suitable to describe both the ground state 
with the ferrimagnetic long rang order and the thermodynamic properties at finite temper- 
atures. The mean field theory gives rise to the leading term in a systematic 1/N expansion 
|| and the effects of fluctuation beyond the mean-field theory can also be discussed, as we 
will mention later. 

In our SBMFT approach, the ground state has a long range ferrimagnetic order arising 
from the condensation of the Schwinger bosons at T=0K. There are two different kinds of 
excitation: One is gapless and ferromagnetic and the other is gapful and antiferromagnetic, 
which has been pointed out by Brehmer et al 0). The spin reduction, the gap of the 
antiferromagnetic branch and the spin correlation at T = OK are calculated and the results 
are in good agreement with those of QMC and DMRG When T > OK, the two 

branches of excitation are both gapful, so that the Bose condensation disappears and there 
is no real long range order. This is just what the Wigner-Mermin Theorem tells us for 
one dimension. The gap of the ferromagnetic branch is proportional to T 2 and the spin 
correlation length is divergent as 1/T. The thermodynamic observables, such as the free 
energy, magnetic susceptibility and the specific heat, are calculated. They agree with the 
numerical results |2],|6j qualitatively and the difference is small at low temperatures. 

This paper is organized as follows: In the second section we construct the mean field 
theory of the ferrimagnetic chain. In the third section we give the ground state properties. 
In the fourth section we study the thermodynamic properties. Conclusions are made and 
advantages and limitations compared with other approaches are discussed in the final section. 

II. SCHWINGER BOSON MEAN FIELD THEORY OF 
THE FERRIMAGNTIC HEISENBERG SPIN CHAIN 

The spin operator Sf can be represented by the Schwinger bosons a^, 
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st = l( a h a i,l -4,1^) , (2) 



with Sf = |(al,T a *,t + a\ ^a^) on each site of kind A. And Sf can be represented in a 
similar way. 

The Hamiltonian(l) is rewritten in this representation 
1 N 

i=l,7/ i,r\ 
N 

= ~ 2 E A,i+?7 A,i+r/ + X! ^i+r/ (^ a ) 
i=l,jj 

Considering the constraint J2a a l,<r a i,<r — 25^ or Y,a^),a b 3,<T — 2S ,S on each site, we may 
introduce two kinds of Lagrangian multipliers \ A and A^ to impose the constraint. At the 
mean field lever, we can take the average value of the bond operator < D iji+ri >= D to be 
uniform and static. And so are < \ A >= \ A and < Af >= \ B . 
The mean field Hamiltonian reads 

N 

H M f = - E { D *( a iA b i+v,l ~ a i,l b i+v,t) + ( a h b l+v,l ~ a tl b i+v,d D } ( 4a ) 

i=l,7j 

Af N 

+A A E K> a ^ - 2SA ) + XB E $,A* - 2SB ) + 2zND*D + zNS A S B , (4b) 

j=l,7J j = l,V 

in momentum space, which is transformed into 
H M f = J2i xAa k,a a k,a + X B bl a b kja } - E{£>^7fcKt & -M - a k,i b -k,i) 

k,a k,a 

+Dz-y k (al A bl ktl - 4,i 6 -m)} + 2zND*D - 2N(S A \ A + S B \ B ) + zNS A S B . (4c) 

Here z is the number of the nearest neighbor sites and equals 2 in the one-dimensional chain, 
and 7fc = - z J2 V e * fcr? — cos k ■ The sum of k is restricted in the reduced first Brillouin zone, 
which extends from — | to |. 

Using the Bogoliubov transformation 
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V 6 -m / 



cosh sinh 
sinh ^fc cosh 



cosh 9k — sinh 
- sinh Ok cosh 6^ 



\ / 3 \ 



with # given by 



tanh 26> 



(A A + A s )/2 



we obtain the energy spectrum 

X A — X B 



X A + \ B . 



\zD~/ k 



E\ 







\ A -\ B 



+ 



,X A + X B . 



2 - \zDj k \ 



k,a 2 ■ y v 2 

The Hamiltonian is diagonalized and it is easy to write down the free energy, 
H MF = £{££(4,^ + 1/2) + El(4,Aa + 1/2)} + 2zND*D 

k,u 

-2N{S A + 1/2)X A - 2N{S B + 1/2)X B + zNS A S B 
F MF 2 n dk ,0 . tn . ,,p 



(5) 
(6a) 

(6b) 
(6c) 



2N (3 



- \l -{ln(2sinh(|^) + ln(2sinh(|Ef)} + zD*D 



A C B 



-(S A + 1/2)X A - (S B + 1/2)X B + ^S A S 



(7a) 



(7b) 



The mean field self-consistent equations can be obtained by minimizing the free energy, i.e. 
SF/SX A = 0, 5F/5X B = and 5F/5D* = 0. After a simple algebra, the equations are 
rearranged as: 



S B -S J 



f -{coth^f-coth^} 
^ 2vr 2 k 2 kf 



s B + S A + 1 = / 



i dk 



\ A +x B 



ft 771/3 1 ft Z7>a1 



2 



{coth + coth '-El} 



z 2n J( * A +* B \2 _ \ z p lk \2 l 2 k 2 " ' 



(8a) 
(8b) 

(8c) 



<A \B 



Rescale the parameters (A , A , D, 0) — > (A 1; A 2 , rj, k) 

A A + A B 1 a n K 

—2— = 2 ZM > = 2 Al7] 

X A — X B 1 n 4k 

= -zA 2 , = — . 

2 2 ' H z 



(9) 



Then the angle of the Bogoliubov transformation is expressed in a compact form 

cosh2# fc = 1 sinh2fl fc = ^ , (10) 



and the self-consistent equations read 
"i dk 



S B -S A = P — {cothKA^l-^ _ A 2 )] _ coth[ K (A lv /l- 77 2 7 2 + A 2 )]} (11a) 
Jon v 



s b + s a +1= n dk coth[ K (A lv /l - ifjj - A 2 )] + coth[/ t (A 1 ^l - gyj + A 2 )] 
7o 7r v/l - ^ 

r- <7£- 

S B + S A + 1 - A^ 2 = / 2 — {cothKA 1A /l-^ 7 2_ A2 )] + 

coth[ K (A 1 ^l-r ? 2 7 | + A 2 )]} v/l - ^ . (11c) 

III.THE PROPERTIES OF THE GROUND STATE 

Notice that the Bogoliubov particles in the (3 branch have to condense to the ground 
state at T = OK as long as S A ^ S B , otherwise the equation (11. a) cannot be satisfied. 
The excitation energy Ep has its minimal value Ep = at k = while the a— branch has 
a gap of 2A 2 at T = OK. Sarker et al [0] showed that the long range order is related to 



the condensation of the Schwinger bosons in both the ferromagnetic and antiferromagnetic 

Heisenberg models. We, now, arrive the same conclusion for the ferrimagnetism model at 

one dimension. This can be contrasted to the antiferromagnetic case, where there is no long 

range order even at T = OK at one dimension. 

Suppose that there is an infinitesimal external stagger magnetic field which is upward at 

the B site and downward at the A site. Then the /3|-branch has the lowest energy and the 

bosons condense at the state of P^k=o- Because of the bose condensation, the self-consist ant 

equations at T = OK are modified as follows: 
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SB ~ SA = 4N coth K Al V 1 - 1 ! 2 - A 2)]U+oo (12a) 



A 1V /1-7? 2 = A 2 (12b) 

S B + S A + 1 = -K( V ) + {SB . ~ (12c) 

S B + S A + 1- A lV 2 = -E(rt) + (S B - S A )Jl-rf (12d) 

7T v 



Here E(rj) and K{rf) are the first and second type complete elliptic integrals. The pa- 
rameters have been determined numerically for the case in S A = 1/2 and S B = 1 with 
7] = 0.8868, Ai = 1.9238, and A 2 = 0.8890. 

The average value of the spin on the site B , 

< Si >= - < 6] T 6 it - b)^ >= — £ < felt & *T - b H b H > , 

k 

can be calculated at T = OK, by using the quasi-particle operator 

< S B Z >= ^-cosh 2 6 k n m \ k=0 ^ 0K = 1(1 + —L=)(S B - S A ) = 0.791. (13a) 

2JM 2 y/1 — I] 2 

Similarly, the average value of the spin on the site A is 

< S A >= -^r- T smh 2 e k n m \ k=0tT -,oK = J(l - J-^ )(S B - S A ) = -0.291. (13b) 

2iV 2 y/1 — T] Z 

The spin reduction r = S B — < S B > on the site B or r = S A + < S A > on the site A is 
given by 

= r* = - V1 _^ 2 l 2 = 0.209 . (13c) 

From equation(lla), we can see that the number of the condensed bosons on the state 
/?! jfc=0 is 2N(S B — S ), which is just the number of the Schwinger Bosons on the site B 
subtracted by that on the site A. As long as S A ^ S B , there is a ferrimagnetic long range 
order, which agrees with Tian's rigorous proof 0. We quote several known results for 
the value of r to show the satisfactory of our calculation. The QMC gives r = 0.207 ± 
0.002, and r = 0.221 in the matrix product states approach 0. The naive SWT overestimates 
the spin reduction and results in r = 0.3 [0. 

The gap of the antiferromagnetic branch in our mean field theory is A ant i = 2A 2 = 1.778, 
which is very close to that in the exact-diagonalization Q A anti = 1.759 and in the QMC 
§, A anti = 1.767. The naive SWT § and the MSWT § give the gap A anti = 1 and 
A anti = 1.676, respectively. 

We calculate the ground state energy of one unit cell, which yields the zero temperature 
free energy per unit cell, 



N 



II £{^1^-^} + \ k \^ ~ 2A ^ SA + S B + l)~ 2A 2 (S A - S B ) + zS A S B 
-1.904 (14) 



This result is much lower than those of the QMC || and the MSWT ||, which are —1.437 
and —1.454, respectively. In fact, this is an artificial result caused by the mean field theory, 
in which we assume that the constraint and the bond operator are uniform and static. We 
have overcounted the degrees of freedom of the Schwinger bosons by a factor of 2, as argued 
by Arovas and Auerbach |||9| . To count the degrees of freedom correctly, we have to divide 
the part of fluctuation per unit cell F MF /N + 2S S B by 2, and add back the classic ground 
energy per cell —2S A S B . Then we have the modified result —1.455. The 1/N expansion 
can give the above argument a strong basis. 

The spin correlation length £ of the ground state is also of interest. Because of the appear- 
ance of the long range order, the transverse and longitudinal fluctuations are anisotrope. In 
our SBMFT, the longitudinal correlation between the site A and the site B can be calculated 



as 



< ^t,z^f,z > < St,z Sf,z > 

= -< (S A Z - at T a iiT )(S£ - b] d b u ) > + < (5& - a\^) >< (S B Z - b) A b jA > 

= - < a h a iAhi > + < a h a iA >< h ),ihi > 

= -\g(R iJ )\ 2 (15a) 
9W) = Jl^ sinh 20 fe e- ifci ^ . (15b) 

Similarly, the longitudinal fluctuations between A site and A site and between B site and B 
site can be given by 

< S&S* >-<S A z >< S A Z >= \f(R tJ )\ 2 (15c) 

< S B Z S B Z >-<S B z >< S B Z >= |/(^)| 2 (15d) 

= — cosh 26 k e- ikR * . (15e) 
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From equations(15.b,15.e) we get the correlation length £ = 77 / [8 ( 1 — 7/ 2 )] 1 / 2 = 0.6785. 
Although the correlation functions are not in a good exponential form because rj is not close 
to 1, we can still see that the correlation decays very rapidly. 
The there kinds of transverse correlation are given by 



< 



-|^-)l 2 - (s B 


- S A )g(R tj . 


) sinh20jfc|fc =o , 


(16a) 


.\f(R ij )\* + 2(S B 


- S A )f{R^ 


)sinh 2 fc | fc=o , 


(16b) 


\f(R ij )\* + 2(S B - 


- S A )f{R i3 ) 


cosh 2 6> fc | fc=0 • 


(16c) 



We note that the transverse correlation length is 2 times of the longitudinal one. The SWT 
calculation gives that the longitudinal correlation length £ = 1/(2 In 2) = 0.7213 0, while 
the numerical methods cannot give the accurate correlation length, because the fluctuation 
decays so rapidly. 

IV. THE THERMODYNAMIC PROPERTIES AT FINITE TEMPERATURES 

We investigate the low temperature asymptotic expansion of the self-consistent equa- 
tions^) and verify that it changes continuously into the equations (12) in T = 0. The 
equation (11. a) has a solution for T 7^ OK. So there is no boson condensation. The equation 
(11. a) gives the gap of the ferromagnetic branch as 

y/l^n 2 2T 2 2 

A f = — ' « 2 4436T 

rf (S B - S^Ax ' 

We solved the self-consistant equations (11) numerically and find that the values of the 
variation parameters are continuously evolved to the values at T=0K. On the other hand, we 
see that the mean field theory fails and the bond operator D is zero when T is higher than 
a specific temperature, say approximately 1.38 in our case. The failure of the mean field 
indicates that the system has entered a local moment phase in which there is no correlation 
of the spin fluctuation, as pointed out by Aromas it et al M. With temperature increasing, 



we find that the gap of the antiferromagnetic branch A anti = A 2 + Ai^l — rf varies. In the 
temperature region where the SBMFT is valid, this gap firstly decreases and then increases. 



It reaches its minimum which is about 1.30 when T rs 0.7. The gaps of the ferromagnetic 
and antiferro magnetic branches are plotted in Figure. 1. And the free energy versa T is 
calculated and plotted in Figure. 2 with the argument of dividing the part of fluctuation by 

2 m- 

The spin correlation at T — > OK are calculated: 



<sf-sf >-- 



1 rf dk 
2J^l 2tt 



2' 



— sinh 2^{coth + coth -E%}e~ ikR >- 



< sf ■ sf >= -m^ 



i.i i 



f dk 



— |...sh 2 9 k coth + sinh 2 9 k coth ^E p k }e~ ikR ^ 



< Sf ■ Sf >= -\F 2 (R 



F 2 (i2y) = f 2 ^{sinh 2 ^ coth S- Et + cosh 2 # fc coth }e 



-i0 \„—ikRi 



(17a) 



(17b) 



(17c) 



Because A/ erro behaves as T 2 , the correlation length behaves as T 1 i.e. £ = rj 2 Ai(S B 

^/(VT^rOT- 1 . 

The dynamic magnetic susceptibility is calculated by using linear response theory, 



K AS^q,u) J 



-g 



*\ <c jS*^, 5^ ^> ^> <c *\ S z S z ^> ^> fl^; 



h A (q,uj) 
h B (q,u) 



(18) 



, where << SfS^ » etc. is the retarded Green function, and h A , h B are the small external 
magnetic field, g is the Laude factor. In the Matsubra representation, the dynamic magnetic 
susceptibilities are given by 



k+q 



-g 2 « SfSf » 



quj 



g* n B (^)-n B (E; 



k+q) 



2N 



iu n + Ei 



El 



cosh 2 6 k+q cosh 2 6^ 



n B {Ej+ q )-n B {El ) _._ , 2n ^ 2n , 1 + n B (E k ) + n B {E\ 
iu n -El +q + Et 



sinh 2 sinh 2 + 



1 (/ cosh 2 sinh 2 9 k 



l + n B (E^)+n B (E p k+q ) 2 . 

COSh 0fc Sinll ; 



fit 3 

n ^k+q 



Ft 



(19a) 



Xab(?,«^) = -si 2 « sfsf » 



2N 



sinh sinh ^fc+q cosh cosh 
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1 + n B {El +q ) + n B {E%) + n B (E% +q ) - n B (E° k 



iu n - Et +q - Et iu n + E% +q - E% 

| n B (El) - n fl (<J 1 + n B (E£ +q ) + n B (Ep 
iu n + El - E p k+q iu n + E% +q + El 

Xbb{Qi i^n) an d Xba(?; i^n) can be obtained by the exchange (a <-> 0) in the equation (19. a) 
and(19.b) respectively. 

The mean field static uniform and staggered magnetic susceptibilities per unit cell are 

..MF r- rib 



N 

.MB 
\-stag 



9 2 P /! + 1) + 4(4 + 1)} (20a) 

9 2 P /! §{[<« + !) + + !)1 ™sh 2 2fl fc + 2 sinh 2 ™k \+"*+J } • (20b) 



If >S A = S B , the above equations are reduced to the familiar forms of the antiferromagnetic 
case |H, which are: 

Auni 2 o 2 uhj 



g 2 (3 f~ 2 —n k (n k + 1) , (20c) 
J=? 2tt 

g 2 (3 / — K(n fc + 1) cosh 2 29 k + k JJ sinh 2 2^} . (20d) 

2iV J=r ^ 13 E 1 ^ 



2N 

Y MF r% rlk 4- 1 /? 

**«er f> ^ {nfeK + l) cosh 2 2^ + ^±i^si 

For low temperatures, both xff n f an d xffog are proportional to T~ 2 , e.g. x«nf = 9 2 A-iV 2 (S B ~ 
S A f/{2^{\ -r/ 2 ))T" 2 , and so on. We plot T Xu ni/Ng 2 versa T in Figure. 3 and T Xs tag/Ng 2 
versa T in Figure. 4. There we have multiply the mean field susceptibilities with the factor | 
due to the same argument as Arovas and Auerbach Txuni/Ng 2 reaches a minimum of 
0.4 at the intermediate temperature region around T »s 0.5. This is due to the contribution 
from the gapful antiferromagntic branch. The low temperature behavior of Txuni/Ng 2 
when T < 0.4 — 0.5, the location and the value of the minimum are found a good agreement 
with with the QMC and DMRG calculations, even better than the MSWT calculation with 
improved dispersion relations . After the point of minimum, the SBMFT result increases 
too rapidly, showing a discrepancy with numerical calculations. Txstag/Ng 2 is dropped 
rapidly and monotonously with the temperature increasing. 

The specific heat C/N versa T is calculated by the numerical differentiation of the 
internal energy with T, and is plotted in Figure. 5. We also find the low temperature 
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behavior of C oc T 1 / 2 , which agrees with the QMC and DMRG calculation |5||| well when 
T < OAK. Again the SBMFT result increases rapidly, failing to see the Shocttky-like peak 
H at intermediate temperatures. 

In short, the SBMFT result describes the thermodynamic properties well at low 
temperatures (T < 0.4 — 0.5). The disagreement in intermediate and high temperatures 
is also owing to the static and uniform constraint. 

V. CONCLUSIONS 

Using the Schwinger boson mean field theory, we have investigated both the ground 
state and the thermodynamic properties of the Heisenberg ferrimagntic spin chain. The 
long range ferrimagnetic order of the ground state is caused by the condensation of the 
Schwinger bosons. The spin fluctuations to the ground state are anisotropic and decays 
very rapidly. The excitation spectrum has both the low energy ferromagnetic branch which 
is gapless at T = OK and the high energy antiferromagnetic branch with a gap is 1.778 
at T = OK. With the temperature increasing, the branch of the ferromagnetic become 
gapful. Af erro behaves as T 2 and A ant i varies. At low temperatures, the ferrimagnetic spin 
chain exhibits the feature of the ferromagnetism. The static magnetic susceptibility, the 
spin correlation length and the specific heat behave as T~ 2 , T~ l and T 1//2 , respectively. At 
intermediate temperatures, the antiferromagnetic branch begins to play an important role 
||. The Txuni has a minimum at T ~ 0.5. 

Compared with other approaches, the SBMFT is a simple mean field theory but can give 
many good results at both zero and finite temperatures. The spin reduction and gap of the 
antiferromagnetic branch A ant i at T = OK differ from the numerical calculations less than 
1%. The thermodynamic properties, such as the static uniform magnetic susceptibility Xuni 
and the specific heat C calculated by the SBMFT, agree with numerical results well when 
T < 0.4. The spin correlation which is anisotropic at the ground state and isotropic at finite 
temperatures can be calculated easily. These results improve those of the SWT largely 
and are consistent with those of the more complicated MSWT and numerical methods ||^| . 
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The SBMFT is not successful at intermediate and high temperatures. The behavior 
of Xuni has quantitative discrepancy with numerical results, and that of C does not agree 
with numerical calculations qualitatively when T > 0.4. Both of them increase too rapidly 
with temperatures increasing. When T > 1.38, the order parameter drops to zero and the 
SBMFT theory fails. It can not describe the system in the whole temperature range as the 
numerical methods and the MSWT do. 
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Figure Caption 



Fig.l The gaps of the ferromagnetic and antiferro magnetic branches 
Fig. 2 The free energy F versa T. 
Fig.3 The Txum/Ng 2 versa T. 
Fig.4 The Txsta g /Ng 2 versa T. 
Fig. 5 The specific heat versa T. 
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Fig. 3 T times the uniform susceptibility per unit site 
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Fig. 4 T times the staggered susceptibility per unit site 
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